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In this paper Einstein's field equations, for static spherically symmetric per- 
fect fluid models with a linear barotropic equation of state, are recast into a 
3-dimensional regular system of ordinary differential equations on a compact 
state space. The system is analyzed qualitatively, using the theory of dynam- 
ical systems, and numerically. It is shown that certain special solutions play 
important roles as building blocks for the solution structure in general. In 
particular, these special solutions determine many of the features exhibited 
by solutions with a regular center and large central pressure. It is also shown 
that the present approach can be applied to more general classes of barotropic 
equations of state. 
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1. INTRODUCTION 

In both Newtonian gravity and general relativity, the simplest models 
of isolated stars are given by static spherically symmetric configurations. 
Despite their simplicity, they are believed to yield many insights about 
much wider classes of stellar models (see, for example, the discussion by 
Hartle [6]). 

The line element for a static spherically symmetric model can be written 

as 

ds 2 = -c 2 ^dt 2 +r(X) 2 [N(X) 2 dX 2 + dn 2 ] , (1) 

with 

dQ 2 = d6 2 + sin 2 6dip 2 , (2) 

where 4>(X) is the gravitational potential, N(\) a dimcnsionless (under 
scale-transformations) freely specifiable function, and r(A) the usual Schwarzschild 
radial parameter, associated with the area of the spherical symmetry sur- 
faces. The coordinate A is a spatial radial variable, whose interpretation 
depends on the choice of N. Since N = 1 corresponds to isotropic coor- 
dinates, the function N can be viewed as a relative gauge function with 
respect to the isotropic gauge. 

This paper is the first in a series devoted to the study of general rela- 
tivistic star models with perfect fluid sources. Thus the energy-momentum 
tensor is assumed to be of the form 

Tab = pU a U b + p (g a b + U a U b ) , (3) 

where p is the energy density, p the pressure, and u a the 4-vclocity of 
the fluid. A relation between the gravitational potential <f) m (1) an d the 
matter content can be found from the equations of motions for the fluid, 
V a T ab = 0, namely 

dp p + p 

This equation is the general relativistic generalization of the corresponding 
Newtonian equation d(j)/dp = — 1/p. 

For a static spherically symmetric perfect fluid model to be considered 
as describing a star, we require that it is isolated in the sense that the 
model has a boundary at a finite radius where the pressure of the fluid 
vanishes. At this radius, the interior solution is matched with a static 
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exterior vacuum spacetime, described by the Schwarzschild solution 

ds 2 = _ (1 _ 2M) rff* + _^!_ + r W , (5 ) 

(see Schwarzschild [12]). A vanishing pressure is, in fact, a necessary and 
sufficient matching condition (see, for example, Stephani [13], p 161). 

In order to analyze the gravitational field equations for (1), an equation 
of state must be specified. Over the years, a number of exact static spher- 
ically symmetric perfect fluid solutions have been found, see, for example, 
Delgaty & Lake [4] for a comprehensive review. Unfortunately, most equa- 
tions of state for known exact solutions have no physical motivation. They 
are instead chosen solely with the purpose of simplifying the differential 
equations, and thereby allowing exact solutions to be found. Our aim is 
not to find new exact solutions, but to gain an understanding of the solu- 
tion space, and its implications, for given equations of state which can be 
physically motivated. 

We will consider barotropic equations of state 

P = p(p) , (6) 

which dominate in the literature (see Stergioulas [14] and references therein). 
Such equations of state are relevant for describing neutron stars and white 
dwarfs, see, for example, Misner et al [2], p 624. As a first step we will con- 
sider complete analytic equations of state covering the pressure range of the 
entire star models. Once the solution space for such equations of states are 
understood, models with composite equations of state (i.e. models where 
the equations of state are obtained by matching different equations of state 
in different pressure regions) can be obtained by matching different matter 
solutions (or in the terminology introduced below, by matching orbits in 
different state spaces). 

One might think that static spherically symmetric perfect fluid models 
lead to simple problems and that there is not much to be discovered in this 
area of research. In this series of papers we will show that this is not the 
case. Some problems turn out to be quite complicated and one can obtain 
new insights. 

The most favored approach for studying static stars, the so-called Tolman- 
Oppenheimer-Volkoff approach, uses the pressure and mass as dependent 
variables and the Schwarzschild radial coordinate as the independent vari- 
able (see, for example, section 6.2 in Wald [18]). The resulting equations, 
however, are not regular at the center of the star, and the problem of prov- 
ing existence and uniqueness of solutions requires a quite technical analysis 
of singular differential equations, see Rendall & Schmidt [11]. 
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Our method is to recast the field equations for a given equation of state as 
a regular 3-dimensional system of ordinary differential equations on a com- 
pact state space. This regularization means that the problems associated 
with the singular nature of the Tolman Oppenheimcr Volkoff equations 
are circumvented. In addition, recasting the field equations into regular- 
ized form on a compact state space allows us to conveniently investigate the 
solution space using powerful methods from the theory of dynamical sys- 
tems. Such methods have been used in spatially homogeneous cosmology 
with great success, see Wainwright & Ellis [17] and references therein. The 
regularized form is also suitable for numerical calculations. The fact that 
the dynamical system is compact and 3-dimensional is of great advantage, 
since this makes it possible to visualize the state space and thereby obtain 
a clear picture of the structure of the entire solution space. 

We will show that there exist special solutions that play important roles 
as building blocks for the remaining solution structure. Some of these 
solutions turn out to determine many of the features exhibited by solutions 
with a regular center and a large central pressure. We will also show that 
the behavior of the solutions at small (and to some extent also large) radii 
is intimately connected with asymptotic self-similarity, even for non-scale- 
invariant equations of state. A useful formulation for the gravitational field 
equations of the static spherically symmetric perfect fluid models and a 
good understanding of the corresponding solution structure may also serve 
as a starting point for exploring their role in a broader context, as will be 
discussed in the concluding remarks. 

In this paper we consider the linear equation of state 

9 = Po + (V ~ 1)P , (7) 

where the constants po and r] satisfy po > 0, 77 > 1. The case r] = 1 corre- 
sponds to an incompressible fluid with constant energy density, while the 
case po = describes a scale- invariant equation of state. The scale- invariant 
case has been investigated previously by Collins [3] using dynamical sys- 
tems theory. Note that w < 2 corresponds to non-causal fluids, in which 
the velocity of sound is greater than that of light. 

The outline of the paper is as follows: in Section 2 we write the grav- 
itational field equations as a 3-dimensional regular dynamical system on 
a compact state space and perform a local analysis of the resulting equa- 
tions. We also present a monotone function and describe the dynamical 
behavior on the boundaries. In Section 3 we focus on the regular solutions, 
and discuss their behavior. In Section 4 we discuss how one can adapt 
the approach of this paper to facilitate studies of more general barotropic 
equations of state. We conclude with some remarks in Section 5. 
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Throughout the paper, geometric units with c = G = 1 are used, where 
c is the speed of light and G the gravitational constant. Roman indices, 
a, b, ... — 0, 1, 2, 3 denote spacetime indices. 



2. THE DYNAMICAL SYSTEM FORMULATION 

In order to recast the gravitational field equations for static spherically 
symmetric prefect fluid models with a linear equation of state into a regular 
dynamical system on a compact state space, we proceed as follows. We first 
write the line element in the form 

ds 2 = -e 2 Ut 2 + d( 2 + e 2 ^ 2 W . (8) 

Then we introduce the variables 

9 = V> , a = 6 , B = , (9) 

where a dot denotes differentiation with respect to £. The gravitational 
field equations, expressed in these variables, are 

6 = -29 2 + 9<r + B 2 + 167rp , (10a) 

& = -20a + a 2 + 4tt(p + 3p) , (10b) 

B = (-0 + a)B , (10c) 

8tt P = 6 2 - a 2 - B 2 . (lOd) 

This system is very similar to those obtained in spatially homogeneous 
cosmology (see Wainwright & Ellis [17]), even though the physical inter- 
pretation is quite different. One can hence import ideas from treatments 
of spatially homogeneous models to the present context. 
To obtain a compact state space we introduce the variables 



{Q,S,C}, (11) 



according to 

Q 



^6 2 + 8TTpo/rj ' s]Q 2 + 87rp /r/ 

B 2 

9 2 + 871750/7/ ' 



R 2 

C = 1^-1— 7- • (12) 



These variables are closely related to those used by Uggla et al [16] in 
cosmology and by Nilsson et al [10] for studying static cylinders with a 
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linear equation of state. We also introduce and a new independent variable 
A, denned by 

- = 1 (13) 

dX ^ + Snpo/ri ■ 

The above choices correspond to 

N2 = c - *' = s - = < 14 > 

in (1), and the prime denotes differentiation with respect to the indepen- 
dent spatial variable A. From (14) it follows that C and 1 — Q 2 are positive. 
Integrating (4), yields 



1-Q 2 ^ 1/v 



l-S 2 -C 



(15) 



where a is a freely specifiable constant corresponding to the freedom of 
scaling the time coordinate t in the line element (1). This, in turn, re- 
flects the freedom in specifying the value of the gravitational potential <fi at 
some particular value of r. Matching an interior solution with the exterior 
Schwarzschild solution, however, fixes this constant. For the purpose of 
interpreting the variables Q and S, it is worth noting the relation 

^ S (16) 



d\nr Q-S 

In terms of the new variables, the gravitational field equations (10a)- 
(10c) takes the form 

Q' = (1 - Q 2 ) (QS - C - 2S 2 ) , (17a) 
5' = \W{2 + r 1 -2QS)-(l-S 2 )(l-Q 2 + QS) , (17b) 
C - 2 [5(1 -Q 2 ) + Q (S 2 -W)]C . (17c) 

Equation (lOd) leads to 

W = 1 - S 2 - C , (18) 

where W is defined by 



W = 8irr]- 1 Cr 2 (p + p) , 



(19) 
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and satisfies W > if we assume that the weak energy condition p + p > 
is satisfied. Equation (18) is used to eliminate W in (17a)-(17c). It follows 
from (18), (17b), and (17c), that 

W = [4QS 2 - 2SQ 2 -r]S + 2QC] W . (20) 

Hence, W — is an invariant subset, as are Q = ±1 and C = 0, which 
is easily seen from (17a) and (17c). These invariant subsets constitute 
the boundary of the physical state space for spherically symmetric models 
with a non-scale-invariant linear equation of state, i.e., they describe the 
boundary of the part of state space where C > 0, W > and po ^ 0. 
We now include these boundaries and obtain a regular dynamical system 
(17a)-(17c) on a compact state space. 

If the variables (11) are used for plane-symmetric models, 1 (17c) de- 
couples and the remaining equations are identical to the C — subset of 
(17a)- (17c). It is therefore natural to refer to this boundary as the plane- 
symmetric boundary. The Q = ±1 subsets correspond to setting p = 
in (7), and the remaining equations describe models with a scale- invariant 
equation of state p = (77 — l)p. 2 We refer to these boundaries as the 
scale-invariant boundaries. The state space, with the different boundaries 
identified, is shown in figure 1. 

An important relation is 

VP _ Q 2 -S 2 -C 

P o~ i-Q 2 ' ( j 

since it yields 

Q 2 - S 2 - C > , (22) 

whenever the pressure is non-negative. The expression Q 2 — S 2 — C = 
defines a surface in state space, which we refer to as the surface of vanishing 
pressure. This surface is also shown in figure 1. It is not an invariant subset 
of (17a)-(17c). 

A monotone function excludes equilibrium points, periodic orbits, recur- 
rent orbits, and homoclinic orbits in its domain. The function 

2Q - S (23) 



V(2Q - sy + 3(1 - Q 2 ) 



lr These models have the same line element as that in (1) but df2 2 = dx 2 + dy 2 . 
2 To obtain the line element for this case, an additional dimensional variable needs to 
be considered as well. 
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which satisfies 

z , = 3(1 -Q 2 ) [2C + V W + 2(2S-Q) 2 ] 
2 [(2Q - S f + 3(1 - Q 2 )] 3/2 

is monotonically decreasing and is defined everywhere in state space ex- 
cept on the scale- invariant boundaries Q = ±1. Hence, the "past" (r — > 0) 
and "future" (r — > +oo) attractors lie on the Q = ±1 boundary sub- 
sets. On these boundaries, however, there exist other monotone functions 
which imply that the only attractors are equilibrium points on the Q = ±1 
boundaries. 3 The equilibrium points of (17a)-(17c), together with their 
corresponding eigenvalues are listed in Table 1. The equilibrium points 
describe self-similar exact solutions to the gravitational field equations. 
The points P^ 2 correspond to the plane-symmetric vacuum solutions, first 
found by Kasner [7]. The points P^ correspond to the Minkowski space- 
time on explicitly spherically symmetric form. The points P± correspond 
to a non-regular self-similar perfect fluid solution due to Tolman [15]. This 
solution is, however, also associated with many other authors, for example, 
Misner & Zapolsky [9]. The points P^ change from focuses to nodes at 
rj = \/64/7 — 2. The equilibrium points Pg^ correspond to a self-similar 
plane-symmetric perfect fluid model. Since this model only exists for r\ < 2, 

3 The equations on the boundary subsets Q = ±1 describe models with a scale- 
invariant equation of state. These equations exhibit a monotone function given in Go- 
liath et al [5]. 
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the fluid is necessarily non-causal. There is thus a bifurcation associated 
with the transition from a non-causal to a causal fluid at 77 = 2. The equi- 
librium points Pc^ leave the state space through the points respectively, 
thereby changing these points from nodes to saddles. 



TABLE 1. 

Equilibrium points and their stability for the linear equation of state 
using the variables {Q,S,C}. The constant a is given by a = 
r) 2 + Ar) — 4. The points P^ 1 only exist for 1 < r\ < 
2. The constant e takes the discrete values ±1. 



Eq point 


Q 


s 


c 


Eigenvalues 




±1 


±1 





±(2 - 77) , ±2 , ±2 




±1 


Tl 





±(6 + 77) , ±6 , ±2 


P ± 


±1 





1 


±2 , T l , ±2 


pi 


±1 


±-2- 

2+J7 




[2 + r7 + eV(2 + -7) 2 -8a] , ±^ 


pi 


±1 


±i(2 + 7?) 





±|(2 + r7) , Ti(2-7?)(6 + 77) , ±±a 



The line element expressed in the variables (11) is invariant under the 
discrete symmetry 

(Q,S,X)^-(Q,S,X) , (25) 

and since (17a)- (17c) are also invariant under this discrete symmetry, a 
solution is represented by two orbits in the state space. We can, however, 
without loss of generality, focus on orbits entering the state space from the 
Q = 1 boundary subset. 

We note that all orbits which correspond to solutions with p ^ 0, with 
spherical as well as plane symmetry, start at equilibrium points on the 
Q = 1 boundary and end at equilibrium points on the Q = — 1 boundary. 
All solutions are thus asymptotically self-similar. However, on their way 
to Q = — 1 from Q — 1 they all intersect the surface of vanishing pressure, 
Q 2 — S 2 — C — 0, at an interior point (Q*, 5*,C») of the state space. To 
obtain physically reasonable spherically symmetric models (models with 
non-negative pressure) one matches each interior solution with the exte- 
rior Schwarzschild vacuum solution at the radius, R, where the pressure 
becomes zero. The radius, R, is determined by inserting (Q*, S 1 *, C*) into 
the expression for r in (14). This leads to 
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From Ql - S% - C* = and C« > 0, it follows that Ql ^ S 2 . Hence all 
spherically symmetric solutions are finite. 

2.1. The boundary structure 

To understand the structure of the interior state space, one has to un- 
derstand the structure of the boundary subsets. Orbits belonging to the 
scale- invariant (Q — 1) boundary are shown in figure 2a, b while orbits be- 
longing to the plane-symmetric (C = 0) boundary are shown in figure 2c, d. 
We refrain from explicitly showing orbits for models with 77 in the interval 
\/64/7 — 2 < 77 < 2, since these models only require a slight modification 
of figure 2a. In addition, the incompressible fluid case rj = 1 (which be- 
longs to the shown 1 < r] < y / 64/7 — 2 interval), and the causal fluid cases 
77 > 2, are the most interesting ones. We require the solutions to have 
non-negative pressure and, as seen from the above discussion, the surface 
of vanishing pressure cut all solutions before they come into the neighbor- 
hood of the W = boundary. Thus this boundary plays a physically less 
important role than the other boundary subsets and we therefore refrain 
from showing orbits belonging to this subset. 



3. REGULAR SOLUTIONS 

Rather than considering all possible orbits, and the associated solutions' 
physical features, we focus on orbits corresponding to solutions with regular 
centers and positive pressure, the so-called regular subset of solutions. 

The spacetime for a regular star is described by the flat Minkowski geom- 
etry at the center. Hence, all orbits belonging to the regular subset start 
from the equilibrium point , which corresponds to Minkowski space- 
time on spherically symmetric form. In the vicinity of P$~ , one finds the 
following approximate expression for the regular orbits 4 

Q = 1 - eie 2A , (27a) 
S = §(e 2 -ei)c 2A , (27b) 

C = I- 4^ 2 c 2A . ( 27c ) 



where t\ and e 2 are small positive constants. It is, however, only the 
quotient < ei/e2 < +00 that parameterizes the 1-parameter subset of 
regular solutions, which can be seen from 

A^-oop Pc (77-lK 2(77 - 1) + (2 + 77) ( £l /e 2 ) ' 1 ' 



4 These expressions are found from the eigenvectors associated with the two positive 
eigenvalues (both equal to 2) of the equilibrium point . 
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(a) 



(b) 





(c) 



(d) 



FIG. 2. Orbits belonging to the boundary subsets for static spherically symmetric 
models with a linear equation of state, using the variables {Q, S,C}, with (a) the scale- 
invariant boundary Q = 1 with 1 < r\ < 1^/64/7 — 2, (b) the scale-invariant boundary 
Q = 1 with r; > 2, (c) the plane-symmetric boundary C = with r\ < 2, and (d) the 
plane-symmetric boundary C = with r; > 2. 



where p c and p c denote the values of the pressure and energy density at 
the center of the star. From (28) we see that there exists solutions with 
a regular center but negative pressure. For solutions with non-negative 
pressure at the center, the quotient ei/e2 is subject to the constraint 



< 



C2 



2 + V ' 



(29) 



where the equality corresponds to a vanishing central pressure. Setting 
€\ = corresponds to an eigendirection in the Q = 1 subset and is associ- 
ated with taking the limit p c — > oo for the interior solutions. The quotient 
Pel pc is a gravitational strength parameter. The Newtonian limit corre- 
sponds to small values of this parameter. Relativistic effects are thus most 
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pronounced when this parameter takes as large values as possible. The 
maximal value ^ = —^r, describing the high pressure limit, is obtained 
when ei = 0, which corresponds to the Q = 1 subset. The Q = 1 subset 
can thus be expected to play a role when probing relativistic effects. 

Solutions describing stars correspond to orbits that start from the equi- 
librium point P^ and satisfy the positive pressure criteria in (29). These 
orbits eventually pass through the surface of vanishing pressure at an in- 
terior point (Q*, S 1 *, C*) of the state space. From (26) it follows that the 
linear equation of state with po ^= leads to star models with finite radii. 
This is expected, since star models always are finite when p ^ as p — > 
(see Rendall & Schmidt [11]). 

The qualitative behavior of the regular subset, projected onto the plane- 
symmetric boundary (C = 0) is shown in figure 3a for 1 < r/ < ^64/7 — 2 
and in figure 3b for models with r\ > 2. The intersection between the orbits 
in the regular subset and the surface of vanishing pressure is indicated by 
the dashed lines. Two orbits constitute the "high pressure" boundary of the 
regular subset. The first orbit belongs to the boundary Q = 1 and starts 
at P3" and ends at . This orbit corresponds to a unique regular scale- 
invariant solution with infinite radius, in contrast to models with po ^ 
(see, for example, Collins [3]). Thus solutions with high central pressure 
are approximately described by this solution near their center. The other 
orbit enters the interior state space from the "Tolman" equilibrium point 
P± . This orbit represents a non-regular solution, but nevertheless plays 
an important role for regular models, particularly those with high central 
energy density and pressure. We will refer to this orbit as the Tolman orbit 
(although a corresponding exact solution is known only for the "stiff" fluid 
case, 77 = 2, see Tolman [15]). 

For the incompressible fluid case and low values of 7], this orbit forms 
the boundary of a "simple" surface formed by the regular subset. For 
larger values of -q (77 > y 7 64/7 — 2), the self-similar Tolman point P^ is 
a focus. This leads to a more complicated situation where "high central 
pressure" orbits in the regular subset spiral around the Tolman orbit, which 
in this case also acts as a "skeleton" orbit for the regular subset. Thus the 
"regular boundary" orbit, the Tolman point P±~, and the Tolman orbit are 
seen to play key roles for the regular subset. Moreover, it is possible to 
obtain approximations for solutions with high central pressure by piecewise 
joining perturbations of the regular boundary orbit and the Tolman orbit. 

In the interior of the star, a mass function m(r) can be defined (see, for 
example, Misner & Sharp [8]). The quotient m(r)/r is dimensionless, and 
can be written as 



m _ C - (Q - Sf 
r 



2C 



(30) 
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(a) 



(b) 



FIG. 3. Orbits belonging to the regular subset for models with a linear equation 
of state, using the variables {Q, S, C}, projected onto the plane-symmetric boundary 
C = for the cases (a) 1 < r\ < 1^/64/7 — 2, and (b) rj > 2. The intersection of the 
regular subset and the surface of vanishing pressure is indicated by the dotted lines. 



Equation (26) and (30) imply that the total mass M of the star is given by 

M = (31) 

This constant is to be identified with the mass parameter of the Schwarzschild 
solution (5) when an interior solution is matched with the exterior Schwarzschild 
solution. 

Besides the orbits of the regular subset and the single orbit starting from 
the equilibrium point P±~, there is also a 2-parameter set of orbits starting 
from the point P^ and a 1-parameter set of orbits starting from P^ (when 
r\ < 2). These two latter sets of orbits, however, describe solutions that 
start out with negative mass. Positive mass is subsequently added and 
the solutions eventually aquire a total positive mass at a sufficiently large 
radius. 



4. COMMENTS ON POSSIBLE GENERALIZATIONS OF 
THE EQUATION STATE 

In this Section we address the issue whether it is possible to modify 
the above formulation so as to be useful for studying models with other 
barotropic equations of state. As will be seen in this series of papers, the 
behavior when p, p — > oo and when p — > constitute key ingredients when 
one attempts to find useful formulations for a given equation of state. The 
latter limit naturally leads to two types of equations of state; those for 
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which p — > when p — > and those for which p — > /?o > when p — > 0. As 
shown by Rcndall & Schmidt [11], the latter case always leads to models 
with finite radii. In the first case the situation is more complicated and 
one sometimes have models with finite radii and sometimes not. 

The present formulation can be modified to cover the situation when 
p — > po > 0, asp — ► 0, and for some equations of state this may lead to a 
useful approach. Let us consider the following class of equations of state: 

P = Po + Hp/po) - l]p , (32) 

where po > and where r\ is an analytic function satisfying 1 < 77 < 00, for 
all non-negative values of the pressure. Thus, 77 is viewed as function and 
not a constant. Equation (21) yields 

(p/po)v(p/po) = Q2 ~f 2 Q 7 C > ( 33 ) 

and hence 

fQ 2 -S 2 -C\ 

To obtain a dynamical system, one first writes the given equation of state 
in the form (32). Then one uses (33) to directly or indirectly determine 
the function 77 in terms of the variables Q, S and C . Finally one introduces 
rj( v i_Q2 ) into the dynamical system (17a)-(17c), possibly changing the 
independent variable in a suitable way. 

This formulation, however, cannot be used when p — > 0,p — > 0, which 
happens for, e.g., polytropic equations of state. In a sequel to this paper we 
will discuss other types of formulations, suitable for treating such problems. 
Moreover, one can also use these formulations when p — > po > 0,p — > 
and treat models with for example, linear equations of state. The present 
formulation has some advantages in this case, however, and it also sheds 
light on these other formulations. It may also be useful to combine the 
various approaches for certain equations of state. 



(34) 



5. CONCLUSION 

We have expressed the field equations as a 3-dimensional regular system 
on a compact state space for static spherically symmetric models with a lin- 
ear equation of state. In this formulation, the existence of regular solutions 
is trivial since they all start from a hyperbolic equilibrium point. The fact 
that all models of this type have finite radii has been naturally incorpo- 
rated into the formalism. We have shown that all models are asymptotically 
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self-similar for small radii. We have obtained a global picture of the solu- 
tion space and this has revealed that certain solutions, the regular scale- 
invariant solution, the self-similar Tolman solution, and the non-regular 
solution associated with the central infinite pressure limit, play key roles 
for understanding the solution structure, and the structure of the regular 
solutions with large central pressure in particular. 

These special solutions exist in all models with linear equations of state, 
including the incompressible fluid case, and basically play the same role, 
although there are some differences depending on if the Tolman equilibrium 
point is a focus or not. For the incompressible fluid, these solutions are 
intimately connected with the Buchdahl inequalities (one of the inequalities 
is an equality for the solution corresponding to the incompressible Tolman 
orbit), see, for example, Buchdahl [1] and Hartlc [6]. In a sequel to this 
paper we will show that similar key solutions exist for other equations of 
state as well. Thus many of the features encountered in this paper are 
typical for large classes of equations of state. 

An interesting application of formulations of this type, for a given equa- 
tion of state, is to probe how physical features depend on, for example, the 
central pressure. In addition, the understanding of the solution space for 
a class of equation of state allows one to investigate how different physical 
features, like, for example, stability properties, depend on the equation of 
state. Another possible application is perturbation theory. Since one has a 
good understanding about the background solutions, one might investigate 
how details of the equation of state affect the possible gravitational wave 
forms 
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